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ABSTRACT 

Aims. We analyze the distribution of the interstellar matter in the environs of the Hn region Gum 31, excited by the open cluster 
NGC 3324, located in the complex Carina region, with the aim of investigating the action of the massive stars on the surrounding 
neutral material. 

Methods. We use neutral hydrogen 21-cm line data, radio continuum images at 0.843, 2.4 and 4.9 GHz, 12 CO(1-0) observations, and 
IRAS and MSX infrared data. 

Results. Adopting a distance of 3 kpc for the H n region and the ionizing cluster, we derived an electron density of 33±3 cm -3 and 
an ionized mass of (3.3±l.l)xl0 3 M Q based on the radio continuum data at 4.9 GHz. The Hi 21-cm line images revealed an Hi shell 
surrounding the Hn region. The Hi structure is 10.0±1.7 pc in radius, has a neutral mass of 1500±500 M Q , and is expanding at 11 
km s _I . The associated molecular gas amounts to (l.l±0.5)xl0 5 M Q , being its volume density of about 350 cnT 3 . This molecular shell 
could represent the remains of the cloud where the young open cluster NGC 3324 was born or could have originated by the shock 
front associated with the H n region. The difference between the ambient density and the electron density of the H n region suggests 
that the H n region is expanding. 

The distributions of the ionized and molecular material, along with that of the emission in the MSX band A, suggest that a photodis- 
sociation region has developed at the interface between the ionized and molecular gas. The copious UV photon flux from the early 
type stars in NGC 3324 keeps the H n region ionized. 

The characteristics of a relatively large number of the IRAS, MSX, and 2MASS point sources projected onto the molecular envelope 
are compatible with protostellar candidates, showing the presence of active star forming regions. Very probably, the expansion of the 
H ii region has triggered stellar formation in the molecular shell. 

Key words. ISM: Hn regions - ISM: individual objects: Gum 31 - stars: early-type - stars: individual: HD 92206 



1. Introduction 

The interstellar medium associated with the birth place of mas- 
sive stars, such as O or early B-type stars, is made up of dense gi- 
ant molecular clouds. Massive stars are characterized by intense 
photon fluxes which ionize and photodissociate the surround- 
ing material. An H n region is a direct consequence of the high 
rate of Lyman continuum luminosity. Initially, the H n region is a 
small, high density region, commonly named ultra-compact H n 
region (UC Hn, Wood & Churchwell 1989). If the photon flux 
rate of the massive star is sufficiently high, the H n region evolves 
into a normal H n region. 

A neutral shell encircles the H n region during the expand- 
ing phase (e.g. Spitzer 1978). The presence of these neutral 
shells is observed in the Hi 21-cm line emission distribution 
(e.g. Deharveng et al. 2003). Molecular line studies have allowed 
the identification of molecular gas following the outer borders of 
H ii regions, indicating the presence of photodissociation regions 
(PDRs). Deharveng et al. (2005) have detected these PDRs in a 
number of H n regions. 

Send offprint requests to: C. Cappa 
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In the present study, we analyze the distribution of the ion- 
ized and neutral material associated with the H n region Gum 3 1 
(Gum 1955) based on H 1 21-cm line emission data, radio contin- 
uum information at different frequencies, and IR and molecular 
data. 

The Hn region Gum 31 is about 15' in size and approxi- 
mately circular in shape (Fig. 1). It is located at (l,b) = (286°12', 
-0°12') in the complex region of Carina and is considered a 
member of the Car OBI association. The SuperCOSMOS image 
(Parker et al. 2005) shows a quite inhomogeneous Hn region, 
with a sharp and bright rim towards lower galactic longitudes 
and lower galactic latitudes, looking fainter and more diffuse to- 
wards higher galactic longitudes. 

Based on data of the radio recombination lines (RRLs) 
H109a and HI 10a at 5 GHz, Caswell & Haynes (1987) found 
that the LSR velocity of the ionized gas in Gum 31 is -18 km 
s , similar to the velocities of other Hn regions in the area of 
the Car OBI association (Georgelin et al. 1986), and derived an 
electron temperature T e — 7100 K. They estimated a flux density 
Ssghz = 35 Jy. 

The excitation sources of the Hn region Gum 31 are the 
OB star members of the open cluster NGC 3324. The bright- 
est star in this cluster is HD 92206, which is the visual double 
star IDS 10336-5806 in the Index catalogue (Jeffers et al. 1963) 
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Table 1. Radio data: relevant parameters. 
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Radio continuum at 4.85 GHz 


Angular resolution 


5:0 


rms noise 


10 mjy beam -1 


Radio continuum at 2.4 GHz 


Angular resolution 


10:4 


rms noise 


12 mjy beam -1 


Radio continuum at 0.843 GHz 


Angular resolution 


43" x 51" 


rms noise 


1 mjy beam -1 


Hi data 


Synthesized beam 


2:4x2:1 


Number of channels 


256 


Velocity coverage 


(-190,+230) km s" 1 


Velocity resolution 


1.64 km s -1 


RMS noise level 


1.6 K 


12 CO(1-0) NANTEN data 


Angular resolution 


2:7 


Velocity coverage 


(-50.+50) km s -1 


Velocity resolution 


0.2 km s -1 


RMS noise level 


1.0 K 



Fig. 1. SuperCOSMOS image showing the Hn region Gum 31. 
The cross marks the position of the multiple system HD 92206. 
The intensity units are arbitrary. 

with a 1 mag fainter companion (HD 92206B) placed 5" to the 
East. Another bright cluster member is located 35" to the SW. 
This star is CD-57°3378, also referred to as HD 92206C in the 
literature. 

The three brightest stars in NGC 3324 have published spec- 
tral types. Both HD 92206A and B are classified as 06.5V, and 
HD 92206C as 09.5V by Mathys (1988). Walborn (1982) classi- 
fies HD 92206A as 06.5V(n) and the component C as 08.5Vp. 

Moffat & Vogt (1975) first carried out photometric observa- 
tions of about 12 stars in the cluster and estimated a color ex- 
cess E(B — V) ~ 0.45+0.05 mag and a distance d - 3.3 kpc. 
Claria (1977), using UBV photometry, confirmed previous re- 
sults by Moffat & Vogt and estimated d — 3.1 kpc. Vazquez & 
Feinstein (1990) derived a distance d — 3.6 kpc for the cluster 
from UBVRI photometry. More recently, Carraro et al. (2001) 
found about 25 new possible cluster members and derived d = 
3.0+0.1 kpc. On the other hand, distance estimates for Car OBI 
are in the range 1.8-2.8 kpc (Walborn 1995). Bearing in mind 
these results we adopted d = 3.0+0.5 kpc for both Gum 31 and 
the ionizing cluster. 

Carraro et al. (2001) find evidence for pre-main sequence 
members beginning at about late B spectral type, which suggests 
an extremely young age for NGC 3324 (< (2-3)xl0 6 yr). The O- 
type members would be stars recently arrived on the Zero Age 
Main Sequence. 

2. Data sets 

2.1. Radio data sets 

We analyzed the radio continuum emission in the region of 
Gum 31 using data obtained at 0.843, 2.42 and 4.85 GHz, 
which were extracted from the Sydney University Molonglo Sky 
Survey (SUMSS) (see Sadler & Hunstead 2001 for a descrip- 
tion), the survey by Duncan et al. (1995), and the Parkes-MIT- 
NRAO (PMN) Southern Radio Survey (see Condon et al. 1993 
for a complete description of this survey), respectively. 



We used H 1 data from the Southern Galactic Plane Survey 
(SGPS) to analyze the neutral gas distribution in the environs of 
Gum 3 1 . These data were obtained with the Australia Telescope 
Compact Array (ATCA) and the Parkes Radiotelescope (short 
spacing information). A Hanning smoothing (Rohlfs 1986) was 
applied to the Hi data to improve the signal to noise ratio. A 
description of this survey can be found in McClure-Griffiths et 
al. (2005). 

The distribution of the molecular material in the region was 
studied using 12 CO(1-0) line data at 1 15 GHz obtained with the 
NANTEN 4 m telescope of Nagoya University at Las Campanas 
Observatory of the Carnegie Institution of Washington, and pub- 
lished by Yonekura et al. (2005). Observational parameters for 
these data were taken from Yonekura et al. (2005). 

The main observational parameters of these databases are 
listed in Table 1 . 

2.2. Infrared data 

We also investigated the dust distribution using high-resolution 
(HIRES) IRAS, and MSX data obtained through /PACfl The IR 
data in the IRAS bands at 60 and 1 00 fim have angular resolutions 
of Ul and U9. The images in the four MSX bands (8.3, 12.1, 
14.7, and 21.3 yum) have an angular resolution of 18'.'3. MSX 
units were converted into MJy ster -1 by multiplying the original 
figures by 7. 133xl0 6 (Eganetal. 1999). 



3. Results 

3.1. Radio continuum emission 

Figure 2 displays the radio continuum image at 843 MHz. The 
image reveals a radio source of 15' in size, coincident in position 
with the H 11 region. There is a remarkable correlation between 
radio and optical emission regions. The strongest radio emission 
region coincides with the brightest section of the optical rim at I 
=* 286° 10'. The region of diffuse emission near (l,b) = (286°20', 



1 IPAC is funded by NASA as part of the IRAS extended mission un- 
der contract to Jet Propulsion Laboratory (JPL) and California Institute 
of Technology (Caltech). 
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Fig. 2. Top: Radio continuum image at 843 MHz. The grayscale 
is from -5 to 180 mJy beam -1 . The contour lines are 5 (5<x), 
20, 50, 100, 150, and 200 mJy beairT 1 . The cross indicates 
the position of HD 92206. Bottom: Overlay of the image at 843 
MHz (contour lines) and the SuperCOSMOS image of Gum 31 
(grayscale). 



-0°4') also has a weak radio counterpart. The fainter optical re- 
gions at (l,b) = (286°13', -0°15') and (286°20', -0°7') also cor- 
relate with weak radio emission. Along with the optical image, 
this radio image shows that the H n region is far from being ho- 
mogeneous and that it is non-uniform in density. Gum 3 1 is also 
detected in the surveys at 4.85 and 2.4 GHz as an isolated radio 
source. The flux densities at 4.85 and 2.42 GHz are 37.7+2.5 
Jy (coincident with the previous estimate by Caswell & Haynes 
[1987]) and 40.7+2.9 Jy, respectively. 

The spectral index of the radio source, as derived from the 
emissions at 2.4 and 4.85 GHz, is a = -0.09+0.20, confirming 
its thermal nature. This result is compatible with the detection of 
RRL at 5 GHz. 

3.2. Hi results 

Circular galactic rotation models predict negative radial veloc- 
ities of up to about -12 km s _1 in the line of sight towards / 



= 286°. However, radial velocities observed in this section of 
the Galaxy are more negative (of up to -30 km s _1 , Brand & 
Blitz 1993), indicating the presence of non-circular motions. 
Consequently, we paid special attention to the neutral gas emis- 
sion distribution at negative radial velocities of up to -50 km 
s" 1 . 

Figure 3 displays a series of H i images within the velocity 
interval from -12.8 to -32.5 km s _1 in steps of 3.3 km s _1 . The 
presence of a region of low H i emission surrounded by regions 
of enhanced emission approximately centered at the position of 
HD 92206 is clearly identified within the velocity range -12.8 to 
-32.5 km s _1 . An almost complete envelope encircles the void. 

The top panel of Fig. 4 displays the H i emission distribution 
within the velocity range -29.3 to -16.5 km s _1 , where the rel- 
atively thick envelope, of about 8'-10' (or 7.0-8.7 pc at d = 3 
kpc), is clearly defined. Within this velocity range, the envelope 
is more easily identified at higher negative velocities, where the 
approaching part of the shell is detected. 

The bottom panel of Fig. 4 displays an overlay of the Hi 
emission distribution and the SuperCOSMOS image of Gum 3 1 . 
The H i envelope clearly anti-correlates with the ionized nebula. 
This envelope is less bright towards higher galactic latitudes. 

The systemic velocity of the structure, which corresponds to 
the velocity at which the feature presents its largest dimensions 
and deepest temperature gradient, is about -23 km s _1 . This ve- 
locity is similar to the velocity of the ionized gas (-18 km s _1 , 
see Sect. 1) as obtained from RRLs, within the errors. 

The morphological agreement between the optical nebula 
and the H i shell and the agreement in velocity between the ion- 
ized and neutral materials indicate that the H i feature is the neu- 
tral counterpart of the ionized nebula. 

3.3. The CO emission distribution 

The distribution of the molecular gas is shown in Fig. 5. The 
upper panel displays the 12 CO(1-0) emission distribution in 
grayscale and contour lines, while the bottom panel displays an 
overlay of the CO contour lines and the optical image of Gum 3 1 . 
The CO gas distribution displayed in the figure was obtained 
from the molecular data cube kindly provided by Y. Yonekura. 
We integrated the 12 CO(1-0) emission within the velocity inter- 
val -27.2 to -14.0 km s _1 . This velocity range is slightly dif- 
ferent from that used by Yonekura et al. (2005) (-30 to -10 km 
s _1 ), since no CO emission was detected for velocities v < -27 
km s _1 , and v > -14 km s _1 associated with Gum 31. Both im- 
ages, that by Yonekura et al. and the one in Fig. 5, are essentially 
the same. 

Intense CO emission regions encircle the brightest sections 
of the optical nebula with CO clumps strikingly bordering the 
bright rim at (286° 10', -0°10'), and near (286°23', -0°15'), 
where the nebula appears diffuse. The CO envelope is open to- 
wards (l,b) = (286°25', -0°5'). The region of relatively faint op- 
tical emission near (l,b) = (286°20', -0°4') coincides with a low 
CO emission region. 

The close morphological agreement between the optical and 
the molecular emissions in the brightest optical region indicates 
that the molecular material is interacting with Gum 3 1 . 

The thick lines in the bottom panel of Fig. 5 delineate the 
C ls O cores found by Yonekura et al. (2005). The dense cores 
coincide with the the brightest 12 CO emission regions. 

The comparison of the molecular gas distribution as shown 
by the 12 CO emission (Fig. 5) with the Hi structure (Fig. 4) 
shows that the neutral envelope around the Hn region has a 
molecular component. Although, both the H i and the molecu- 
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Fig. 3. Series of Hi brightness temperature images for the velocity interval -12.8 to -32.5 km s _1 averaged over 3.3 km s . The 
greyscale is from 70 to 120 K for the images at v > -26.0 km s _1 , and from 10 to 60 K for the images at v < -26.0 km s -1 . The 
contour lines are from 10 to 1 10 in steps of 10 K. The star indicates the position of HD 92206. The velocity interval corresponding 
to each image is indicated. 



lar envelopes are approximately coincident, most of the bright 
molecular clumps anticorrelate with the H i maxima, indicating 
that gas in these regions is mainly in molecular form. However, 
the comparison of the CO distribution with the H i emission dis- 
tribution at different velocities shows that the 12 CO clump at (l,b) 
= (286°23', -0°15') coincides with an H i cloud present in the ve- 
locity interval from -26 to -29.3 km s _1 . A possible explanation 
is the presence of the CO clump inside a dense H i cloud. 

The comparison between the ionized, neutral atomic, and 
molecular distributions around Gum 3 1 suggests the presence of 
a prominent PDR bordering the brightest ionized regions. 

3.4. The emission in the infrared 

The distribution of the IR emission at 60 and 100 /mi, due to 
thermal dust emission, is displayed in Fig. 6. The upper panels 
show the emission at both wavelengths while the bottom pan- 
els display overlays of the IR emission and the optical image 
of Gum 3 1 . The images reveal an IR structure which is brighter 
near (l,b) = (286°23',-0° 15') and where the optical nebula has its 
sharpest border, and is weaker on the opposite side. The bright- 
est IR emission regions at 60 //m and 100 /mi are projected 
onto the neutral envelope, delineating the ionized nebula. The 
IR emission at (286°23',-0°15') detected at both wavelengths 



coincides with strong 12 CO(1-0) emission. 12 CO emission also 
appears bordering the two IR clumps located near the bright rim 
at (l,b) = (286°10', -0°10') detected at 60 /an, and in between. 
The distribution of the IR emission at both wavelengths shows 
the presence of dust most probably related to the surrounding H i 
and molecular shells. 

Following the procedure described by Cichowolski et al. 
(2001), we derived the color temperature of the dust associated 
with the H n region and the neutral envelope based on the IR 
fluxes at 60 /mi and 100 /mi. Taking into account different val- 
ues for the background emission, we found Td = 34+7 K. The 
range of temperatures corresponds to n = 1-2 and to different IR 
background emissions. The parameter n is related to the dust ab- 
sorption efficiency (k v oc y"). We adopted k v — 40 cm 2 g -1 . This 
value was derived from the expressions by Hildebrand (1983) for 
n = 1 and A = 60 /mi. The dust temperature obtained is typical 
for H ii regions. 

Figure 7 shows an overlay of the distribution of the emission 
in the MSX bands A and E, and the optical image. The emission 
in band A closely follows the brightest sections of the nebula, 
and correlates with the neutral gas. Particularly, the brightest re- 
gions emitting at 8.3 /mi coincide with the dense cores 2 and 6 
(see Fig. 5) found by Yonekura et al. (2005). On the contrary, 
the emission in band E appears clearly associated with the ion- 
ized gas. Note that the strongest emission region in band A at 
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Fig. 6. Top: Far infrared IRAS images at 60 and 100 pm towards Gum 31. The grayscale is from 300 to 3000 MJy ster~' for the 
image at 60 pm, and from 600 to 4500 MJy ster~' for the image at 100 /im. The contour lines are from 200 to 1000 MJy ster~' 
in steps of 200 MJy ster~' and from 1000 to 3000 MJy ster -1 in steps of 500 MJy ster -1 for the image at 60 pm, and from 400 to 
1000 MJy ster -1 in steps of 200 MJy ster~' and from 1000 to 3000 MJy ster~' in steps of 500 MJy ster~' for the image at 100 pm. 
Bottom: Overlay of the optical (SuperCOSMOS) (grayscale) and the IR images. 



(l,b) = (286°23' ,-0°15') coincides with bright molecular and far 
infrared clumps. 

The emission distribution in band A is most probably related 
to emission from polycyclic aromatic hydrocarbons (PAHs). 
According to Cesarsky et al. (1996), these dust grains cannot 
survive inside the H n region, but can on the neutral PDR, where 
they radiate in the PAH bands at 7.7 and 8.6 /im, included in the 
MSX band A. MSX band E, on the contrary, includes continuum 
emission from very small grains, which can survive inside ion- 
ized regions (cf. Deharveng et al. 2005), and a contribution from 
nebular emission lines. 

In summary, the distribution of the neutral atomic and molec- 
ular gas, and that of the interstellar dust reveals the presence of 
a neutral shell surrounding the H n region. The distribution of 
the molecular gas and that of the emission in the MSX band A 
related to PAHs strongly suggests the presence of a PDR at the 
interface between the ionized and the molecular gas. 

Both the infrared emission at 60 and 100 pm and the optical 
image suggest that the faint optical emission region near (l,b) = 
(286°20', -0°4') (indicated by an arrow in the bottom right im- 
age of Fig. 6) are also linked to Gum 3 1 . As pointed out in Sect. 
3.1, this region is also faint in the radio continuum (see Fig. 2), 
and corresponds to weak regions in the H i and CO envelopes 
suggesting that the ionizing flux of the massive stars in the open 
cluster may drain through these regions to the general interstel- 
lar medium. The situation resembles the case of the stellar wind 
bubble around WR 23 (Cappa et al. 2005). 



4. Discussion 

4. 1 . Physical parameters 

The main parameters of the dust and the ionized and neutral gas 
related to Gum 3 1 are summarized in Table 2. The parameters 
of the ionized gas were derived from the image at 4.85 GHz. 
The uncertainty in the flux density corresponds to an error of 0.1 
Jy beam -1 in the estimate of the radio continuum background. 
The electron density and the H n mass were obtained from the 
expressions by Mezger & Henderson (1967) for a spherical Hn 
region of constant electron density (rms electron density n e ). The 
presence of He n was considered by multiplying the H i mass by 
1.27. The number of UV photons necessary to ionize the gas 
log NLy- c was derived from the radio continuum results. Errors in 
the linear radius, in the rms electron density, and in the excitation 
parameter come from the distance uncertainty. The high electron 
density of the H n region is compatible with the relatively short 
lifetime for the cluster. 

The parameters of the neutral atomic gas includes: the (l,b) 
position of the centroid of the Hi shell, the velocity interval 
spanned by the structure, the systemic and expansion velocities, 
the radius of the neutral gas structure, and the associated atomic 
mass. 

The expansion velocity was estimated as in previous papers 
(see Cappa et al. 2005 and references therein) as v exp = (v2-Vi)/2 
+ 1 .6 km s _1 . The extra 1 .6 km s _1 allows for the presence of H i 
in the caps of the expanding shell, which should be present in 
the images at vi + 1.6 km/s and V2-I.6 km/s. These caps are not 
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Fig. 7. Overlay of the MSX infrared images (contours) corresponding to bands A (8.28 yum) and E (21.34 fim), and the 
SuperCOSMOS image of the nebula (grayscale). The contour lines are 25, 39, 57, 85, 1 14, and 140 MJy ster 1 for band A; and 36, 
46, 57, 85, 1 14, and 140 MJy ster 1 for band E. 



Table 2. Main parameters of the ionized and neutral gas in Gum 3 1 



Adopted distance (kpc) 


3.0+0.5 


Hn region 




Flux density at 4.85 GHz (Jy) 


37.7+2.5 


Angular radius (arcmin) 


7.5±0.2 


Linear radius (pc) 


6.5+1.0 


Emission measure (pc cm -6 ) 


(1.5±0.2)xl0 4 


rms electron density [;j e ](cirr 3 ) 


33±3 


Used Lyman UV photons [logN Ly _ c ] (s~') 


49.0±0.1 


Ionized mass (M Q ) 


3300±1100 


Neutral atomic shell 




(l,b) centroid 


286° 15' -0° 10' 


Velocity interval vi,V2 (km s~') 


-13,-32 


Hi systemic velocity (km s~') 


-23 


Expansion velocity v exp (km s~') 


11 


Radius of the H i structure (arcmin) 


11.5 


Radius of the H i structure R (pc) 


10.0±1.7 


Atomic mass in the shell (M Q ) 


1500+500 


Molecular shell 




Velocity interval vi,V2 (km s~') 


-27.2,-14.0 


Mean H 2 column density (cm -2 ) 


1.2X10 22 


Molecular mass (M Q ) 


(l.l±0.5)xl0 5 


Dust related to the Hn region and the neutral shell 




Total dust mass (M Q ) 


60±20 


Dust color temperature (K) 


34±7 



detected in the present case, as in most of the cases, probably 
because of their weak brightness temperature. 

The radius of the Hi structure was estimated from Fig. 4 
and corresponds to the position of the maxima in the shell. The 
neutral atomic mass corresponds to the mass excess in the shell, 
assuming that the gas is optically thin and including a He abun- 
dance of 10%. 

The Hi column density (Nm) and the molecular mass were 
estimated from the 12 CO data, making use of the empirical rela- 
tion between the integrated emission Wco (= jTdv) and Nhi- 
We adopted N H2 = (1.06±0.14) xW C0 x 10 2() crrT 2 (K km 
s -1 ) -1 , obtained from y-ray studies of molecular clouds in the 
Orion region (Digel et al. 1995). To derive the molecular mass, 
we integrated the CO emission within a circle of about 16' in 
radius (= 15 pc at 3.0 kpc), centered at (l,b) = (286°16',-0°8'). 



The ambient density derived by distributing the molecular 
mass over a sphere of 15.0 pc in radius is ^350 crrT 3 , reinforcing 
the idea that the molecular gas represents the remains of the orig- 
inal material where the open cluster NGC 3324 was born. The 
difference between the electron density and the ambient density 
indicates that the H n region is expanding. 

It is known that massive stars have stellar winds with high 
mass loss rates and terminal velocities (Chlebovski & Garmany 
1991; Prinja et al. 1990), which generally affect the surrounding 
interstellar medium creating interstellar bubbles (Weaver et al. 
1977). In the optical domain, they appear as ring-shaped struc- 
tures (Lozinskaya 1982). The optical appearance of the Hn re- 
gion, almost spherical without clear evidence for a central cavity, 
suggests that the massive stars in the cluster have weak stellar 
winds (which is supported by the short age of the cluster) and/or 
have existed during a very short period of time to create an in- 
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Fig. 4. Top: H i brightness temperature image corresponding to 
the velocity interval -29.3 to -16.5 km s _1 showing the H i struc- 
ture related to Gum 31. The grayscale is from 40 to 100 K. The 
contour lines are from 40 to 100 K in steps of 10 K. Bottom: 
Overlay of the SuperCOSMOS (grayscale) and Hi (contours) 
images showing the close correspondence of neutral and ionized 
gas emissions. 



terstellar bubble in an interstellar medium as dense as observed 
here. 

Figure 8 shows schematically the distribution of the different 
gas components. Although Hi and molecular clumps anticorre- 
late in position, the large scale H i and molecular gas coincide 
in position around the nebula but in the region at higher galactic 
latitudes, where the CO emission is observed outside the H i re- 
gion. The large amount of molecular gas compared to the neutral 
atomic gas supports the idea that the CO emission which encom- 
passes the ionized nebula represents the remains of the molecu- 
lar material where the open cluster was born. Part of this H i gas 
may have originated in the photodissociation of the molecular 
gas. 

4.2. Energy budget 

Bearing in mind that massive stars have a copious UV flux ca- 
pable of ionizing the surrounding H i gas, we investigate in this 
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Fig. 5. Top: 12 CO(1-0) emission distribution corresponding to 
the velocity interval -27.2 to -14.0 km s _1 . The grayscale 
is from 15 to 150 K km s _1 . The contour lines are 20.2, 
40.5, 60.7, 81.0 and 101.2 K km s -1 . Bottom: Overlay of the 
SuperCOSMOS (grayscale) and the CO (contours) images. The 
thick lines delineate the C ls O cores described by Yonekura et 
al. (2005). The numbers correspond to the C 18 cores as identi- 
fied in Table 2 by Yonekura et al., and the crosses mark the core 
positions indicated in the same table. 



section whether the massive stars in NGC 3324 can provide the 
energy to ionize the gas. 

As described in Sect. 1, the only O-type stars in the open 
cluster are the multiple system HD 92206, which contains three 
O-type stars classified as 06.5V, 06.5V and 08.5 (Walborn 
1982, Mathys 1988). Considering that HD92206A is almost 1 
mag brighter than component B (Claria 1977), the similar spec- 
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Table 3. IR Point sources and YSO candidates from the IRAS, MSX and 2MASS catalogs 



# 


l[o '] 


b[o '] 


IRAS source 




Fluxes [Jy] 




L iras [WLq] 










12/jm 


25//m 


60//m 


100/jm 




1 


285°59:82 


+0°05:82 


10349-5801 


1.6 


2.5 


27.9 


104 


1.5 


2 


286°04:98 


-0°25:08 


10335-5830 


0.8 


3.5 


22.8 


111 


1.5 


3 


286°08:58 


-0°06:96 


10351-5816 


6.3 


6.8 


315 


1480 


18 


4 


286°08:64 


-0° 18:78 


10343-5826 


5.5 


4.3 


103 


340 


5 


5 


286° 10:92 


-0°14:7 


10349-5824 


5.9 


9.0 


163 


1660 


18 


6 


286 12: 18 


+0 10.20 


10365-5803 


7.2 


86 


1170 


2780 


43 


/ 


ZOO lj.12 


n°oc A A 
— U Zj .44 


1 HQ A £. COTC 

1U34d-jo3j 


t i 
1.1 


3.3 


11.7 


1 A Q A 

143U 


1 A 

14 


8 


286 17:34 


+0 00:42 


10365-5814 


2.3 


1.7 


85 


246 


4 


9 


286 22:5 


-0°15:3 


10361-5830 


12.4 


38.4 


626 


2160 


30 


10 


286 26:94 


-0 22:92 


10361-5839 


2.8 


5.4 


51.6 


307 


4 


11 


286 29:34 


-0° 16:68 


10368-5835 


2.1 


2.3 


84 


271 


4 


12 


286 33:72 


ao Any An 

-0 07:08 


10379-5828 


1.0 


2.0 


13.8 


540 


5 


# 


/[C] 


b[o '] 


MSX source 




Fluxes [Jy] 














8/jm 


12//m 


14yum 


21/jm 


Class. 


13 


286°09:78 


-o° 11:28 


G286. 1626-00. 1877 


0.7311 


1.298 


1.256 


2.806 


CHu 


14 


286° 12!36 


-0°09:66 


G286.2056-00.1611 


0.1585 


0.9313 


2.151 


6.904 


MYSO 


i c 


2S6 12.4a 


no 1 
— U 11). 32 


pig/: AA 1 -7TA 

O2S0.2U / /-UU. 1/21) 


A AO-7A 
U.U8 /U 


0.9481 


1.764 


O A 

2.y 


nvcn 
MYi>(J 


16 


286° 12:54 


+o°io:i4 


G286. 2086+00. 1 694 


1.353 


2.882 


7.182 


40.57 


MYSO 


17 


286° 12:6 


-0° 10:68 


G286. 2096-00. 1775 


0.2202 


0.9538 


1.368 


7.206 


MYSO 


18 


286° 15:42 


-0° 19:44 


G286. 2566-00.3236 


2.04 


2.151 


1.34 


4.328 


CHu 


19 


286°2i:48 


-0° 17:58 


G286. 3579-00.2933 


0.7126 


1.815 


2.677 


6.065 


MYSO 


20 


286 22:5 


-0° 15:78 


G286. 3747-00.2630 


3.591 


4.756 


2.409 


7.577 


CHu 


21 


286°22:62 


-0° 15:36 


G286. 3773-00.2563 


1.628 


2.918 


3.855 


12.1 


MYSO 


# 


l[o'] 


b[o>] 


2MASS source 


-/[mag] 




/(' lm q it 1 


(J-H) 


(H-K) 


22 


285°54:42 


-o° 13:74 


10350210-5831039 


10.584 


lU.Dj 1 


1U.JCH- 


-0.047 


0.067 


23 


286°04:44 


-0=00:24 


10365972-5824186 


12.042 




y.zoH- 


1.532 


1.226 


24 


286°06:36 


-0°08:64 


10364112-5832326 


11.367 


11.385 


11.28 


-0.018 


0.105 


25 


286°06:54 


-0°08:40 


10364296-5832267 


10.465 


10.513 


10.439 


-0.048 


0.074 


26 


286°09:60 


-o°oo:30 


10373406-5826540 


12.487 


11.454 


10.513 


1.033 


0.941 


27 


286°09:78 


-o° 11:22 


10365396-5836293 


12.242 


10.99 


10.103 


1.252 


0.887 


28 


286°10:08 


-0°09:i2 


10370395-5834489 


10.772 


9.712 


8.685 


1.06 


1.027 


29 


286° 10:20 


-o°n:io 


10365749-5836366 


13.99 


12.815 


11.984 


1.175 


0.831 


30 


286° li: 16 


-0°02:64 


10373574-5829405 


15.318 


13.343 


11.614 


1.975 


1.729 


31 


286° 12:96 


-0°04:86 


10373956-5832311 


10.926 


10.752 


10.554 


0.174 


0.198 


32 


286° 13:02 


-o° 10:86 


10371717-5837460 


11.793 


11.773 


11.665 


0.02 


0.108 


33 


286° 1332 


-0°08:46 


10372824-5835492 


12.191 


11.301 


10.448 


0.89 


0.853 


34 


286° 13:38 


-o° 10:20 


10372226-5837229 


7.563 


7.588 


7.479 


-0.025 


0.109 


35 


286° 13:92 


-0° 17:87 


10365763-5844052 


12.22 


11.882 


11.51 


0.338 


0.372 


36 


286° 19:26 


-0°06:84 


10381421-5837192 


12.635 


11.886 


11.349 


0.749 


0.537 


37 


286°20:i0 


-0° 19:80 


10373105-5849026 


15.155 


12.8 


11.127 


2.355 


1.673 


38 


286°2i:06 


-0° 16:86 


10375219-5847133 


12.271 


11.363 


10.675 


0.908 


0.688 


39 


286°2i:90 


-0°04:92 


10383875-5836566 


12.116 


11.636 


11.187 


0.48 


0.449 


4U 


2S6 22.92 


no 1 Q /OA 

— U 13. 21) 


11)38 140 1-5844410 


12.36 / 


11.815 


11.398 


ft rcT 

U.332 


A A 1 H 

(J. 41 / 


41 


286° 24: 92 


-0° 18:66 


10380736-5850240 


12.477 


11.538 


10.78 


0.939 


0.758 


42 


286°28:il 


-0°24:50 


10380702-5857039 


12.724 


11.599 


10.726 


1.125 


0.873 


43 


286°28:44 


-0°23:71 


10381226-5856318 


13.348 


12.274 


11.355 


1.074 


0.919 


44 


286°29:46 


-0°06:30 


10392451-5841486 


13.836 


12.561 


11.619 


1.275 


0.942 


45 


286°29:46 


-o°08:io 


10391799-5843257 


13.195 


12.322 


11.624 


0.873 


0.698 


46 


286°30:90 


-0°06:30 


10393410-5842321 


11.701 


10.423 


9.334 


1.278 


1.089 


47 


286°3i:27 


-0°28:40 


10381363-5902003 


13.368 


12.557 


11.982 


0.811 


0.575 



tral types found by Mathys (1988) suggest that HD92206A is 
probably a spectroscopic binary 06.5V + 06.5V. Taking into 
account the UV photon fluxes emitted by the stars, N* (s _1 ), de- 
rived by Martins et al. (2005) from stellar atmosphere models, a 
group of four massive stars having the spectral types indicated 
above have a total UV photon flux corresponding to log N* = 
49.4. By comparing N* with the UV photons used to ionize the 
gas ( NLy-c) listed in Table 2, we conclude that the massive O- 
type stars in NGC 3324 can maintain the H n region ionized. 

The radius of the Stromgren sphere formed in a region with 
a volumetric ambient density of 350 cirT 3 is ^ 1.7 pc, which is 
lower than the radius of the ionized region. Following Spitzer 



(1978), we estimated that the H n region has been expanding for 
*lxl0 6 yr. 

The ambient density we derived from the 12 CO(1-0) data is 
larger than the rms electron density, 33 crrT 3 (see Table 2), giv- 
ing additional support for the interpretation that the H n region 
is expanding. 

5. Stellar formation 

Stellar formation may be induced by the expansion of the H n 
region in the surrounding molecular envelope, where the pres- 
ence of high density regions favors the conditions which lead 
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GALACTIC LONG. 

Fig. 9. Top: Point sources from the IRAS (triangles), MSX (circles), and 2MASS catalogs overlaid onto the SuperCOSMOS image 
and the 12 CO and C ls O contours (see Fig. 5 for details). 

to the formation of new stars. Among the different processes 
that induce star formation, the collect and collapse process pro- 
posed by Elmegreen & Lada (1977) may work efficiently in the 
dense molecular shells around H n regions. The physical condi- 
tions in the molecular shells were discussed by different authors 
(Whitworfh et al. 1994; Hosokawa & Inutsuka 2005; 2006, Dale 
et al. 2007). 

To investigate the presence of protostellar candidates in this 
region we used data from the IRAS, MSX and 2MASS point 
source catalogs. We searched for point sources in a region of 
about 20' in radius centered at the position of the NGC 3324. 

5. 1 . IRAS sources 

The IRAS point source catalog allows identification of proto- 
stellar candidates following the criteria by Junkes et al. (1992). 
Sources with quality factors (260 + 2ioo ^ 4 were considered. 
Twelve out of the thirteen point sources have IR spectra com- 
patible with protostellar objects. The names of the IRAS proto- 
stellar candidates, their (l,b) position, their fluxes at different IR 
wavelengths, and their luminosity derived following Yamaguchi 
et al. (2001), along with a reference number are listed in Table 
3. The IRAS protostellar candidates are indicated as filled trian- 



gles in Fig. 9, superimposed onto the molecular gas and ionized 
gas distributions. Each source can be identified in the figure by 
its reference number. They are probably young stellar objects 
(YSOs). 



5.2. MSX sources 

Massive young stellar objects (MYSOs) can be identified from 
the MSX point catalogue following the criteria by Lumsden et 
al. (2002). A total of 310 MSX point sources were found to be 
projected onto the area. 

Lumsden et al. (2002) found that MYSOs have infrared 
fluxes with ratios F2\/F% > 2 and Fu/Fu > 1, where Fg, Fu, 
F14, and F21 are the fluxes at 8.3, 12, 14, and 21 p.m. For com- 
pact Hn regions, ratios are Fi\IF% > 2 and F^/F^ < 1- Taking 
into account sources with flux quality q > 2, we were left with 14 
sources after applying Lumsden et al.'s criteria indicated above. 

Six out of the 14 sources can be classified as MYSOs, while 
three sources are compact H 11 regions, also indicative of active 
stellar formation. The nine sources are listed in Table 3 and in- 
dicated in Fig. 9 as circles. 
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Fig. 8. Schematic view of Gum 3 1 and the respective locations 
of the different gas components and of interstellar dust. 



5.3. 2MASS sources 

Point sources with infrared excess, which are YSO candidates, 
were searched for in the 2MASS catalog (Cutri et al. 2003), 
which provides detections in three near IR bands: J, H, and 
K s , at 1.25, 1.65, and 2.17 fim, respectively. A total of 2xl0 4 
sources are projected onto a circular region of 20' in radius. 
We took into account sources with S /N > 10 (corresponding 
to quality "AAA"). Only sources with Ks < 12 were included. 
The last criterium corresponds to stars with spectral types ear- 
lier than B3 at a distance of 3.0 kpc. Following Comeron et al. 
(2005) and Romero (2006), we determined the parameter q = 
(J -H)- 1.83 x (H - K s ). Sources with q < -0.15 are classi- 
fied as objects with infrared excess which may be YSOs. After 
applying this criterium we were left with 26 sources, which are 
shown in the color-color and magnitude-color diagrams of Fig. 
10. Absolute magnitudes for the ZAMS in the Ks band are from 
Hanson et al. (1997), while values corresponding to (J-H) and 
(H-K s colors were obtained from Koornneef (1983). The loca- 
tion of these sources is also marked in Fig. 9 as crosses. The 
main data of these sources are shown in Table 3. 

The selected sources are located to the right of the redden- 
ing vector in the color-color diagram, where YSOs are expected 
to be placed. Most of the sources with the highest infrared ex- 
cess appear projected onto the molecular clouds (see Fig. 9). 
The magnitude-color diagram shows that most of the sources are 
MYSOs. Although these diagrams are not conclusive in identify- 
ing YSOs, the strong infrared excess of the sources is compatible 
with protostellar candidates. 

5.4. Distribution and characteristics of the YSOs 

Figure 9 shows that the IRAS and MSX point sources, and most 
of the 2MASS point sources classified as YSOs appear border- 
ing the ionized region, projected onto the molecular envelope 
detected in 12 CO emission, close to the periphery of the Hn re- 
gion, or near the central cluster. Some are also coincident with 



the dense cores found by Yonekura et al. (2005) in C O. We 
analyze some particular regions. 

The IRAS source #9, the MSX sources #19, #20, and #21, 
and the 2MASS sources #38 and #40, are projected onto a 12 CO 
clump and onto the dense core 6 found in C 18 0. The presence of 
these sources indicates that stellar formation is ongoing in this 
particular molecular clump. As suggested by Yonekura et al., a 
star cluster including massive stars is probably being formed in 
this region. The MSX source #19 is coincident with an extended 
radio source at 843 MHz (see Fig. 2) and with a bright source in 
the MSX band E (see Fig. 7), showing that a compact H n region 
may be present inside the molecular cloud. The MSX sources 
#20 and #21 coincide with a small source detected in both bands 
A and E (see Fig. 7) and with a point radio source detected at 
843 MHz (see Fig. 2). The detection of radio continuum emis- 
sion is compatible with the classification of compact H n region 
of source #20. The fact that IR emission was detected in both 
bands A and E suggests contributions of warm interstellar dust 
and/or ionized gas. The radio source coincident with source #19 
is about 48"x38" in size and its flux density is about 30 mJy. No 
electron density was estimated since, very probably, the compact 
H ii region is optically thick at this low frequency. These results 
reinforce the idea that this is a very active star forming region. 

Particularly interesting is the bright rim region at (l,b) - 
(286°10',-0°10'). The 2MASS sources #27, #28, and #29, and 
the MSX source #13 are located in this area, projected onto the 
molecular envelope. IRAS sources #3 and #5 are also placed 
close to the border of the H n region, coincident with molecu- 
lar material. Source #3 coincides with a region emitting in the 
MSX band A and with core 2. 

A bunch of protostellar objects is almost coincident in po- 
sition with the open cluster NGC 3324: three MYSOs (MSX 
sources #14, #15, and #17), and the 2MASS sources #32, #33, 
and #34. The color-color diagram shows that two of these 
sources have relatively small IR excess. Some of these sources 
coincide with the loose IR cluster IC2599 listed by Dutra et 
al.(2003). These facts suggest that stellar formation is ongoing 
in the region of NGC 3324, as previously found by Carraro et 
al. (2001). Also IRAS sources #6 and #10 are projected onto the 
dense cores 5 and 7, respectively. The large luminosities L IRAS 
estimated for some of these sources suggest that they are proto- 
stellar candidates for massive stars or star clusters. 

The sources IRAS 10355-5828 (286° 16:92,-0° 15:6) and 
MSX G286.2868-00.2604 (neither included in Table 3 nor 
shown in Fig. 9) coincide with HD 92207, a red supergiant 
whose membership to the open cluster is a matter of debate 
(Carraro et al. 2001). Its MSX fluxes correspond to an evolved 
object. 

The detection of protostellar candidates in the IRAS, MSX, 
and 2MASS databases, strongly indicates that active stellar 
formation is currently ongoing in the molecular shell around 
Gum31. 

To sum up, most of the protostellar candidates detected to- 
wards Gum 3 1 appear projected onto the shell detected in the 
12 CO line, and coincide with the dense cores detected in C lo O 
emission. Some are located close to the bright sharp borders of 
the H ii region, and near the open cluster NGC 3324. 

The presence of protostellar objects on the molecular enve- 
lope bordering the ionized region indicates that star formation 
has been triggered by the expansion of the H n region. The dis- 
tribution of the molecular and H i gas around the ionized region 
suggests that star formation could be due to the collect and col- 
lapse process. 
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Fig. 10. Color-color and magnitude-color diagrams showing the 2MASS sources with infrared excess. 



6. Summary 

We have analyzed the interstellar medium in the environs of the 
H ii region Gum 3 1 to investigate the action of the massive stars 
in the exciting open cluster NGC 3324 on the surrounding neu- 
tral material. 

We based our study on H i 21 -cm line emission data belong- 
ing to the SGPS, radio continuum data at 0.843, 2.4 and 4.85 
GHz from the PMN Southern Radio Survey, I2 CO data from 
Yonekura et al. (2005), and IRAS (HIRES) and MSX data. 

Adopting a distance of 3.0+0.5 kpc, we have derived an ion- 
ized mass of 3300±1100 M and an electron density of 33+3 
cm 4 . The four O-type stars in the HD 92206 multiple system 
can provide the necessary UV photon flux to maintain the H n 
region ionized. 

The Hi emission distribution in the environs of Gum 31 
shows the presence of an H i shell approximately centered at the 
position of the multiple system HD 92206. The H i shell closely 
encircles the optical nebula. It is detected within the velocity 
range -32 to -13 km s -1 and its systemic velocity of -23 km s _1 
is coincident, within errors, with the velocity of the ionized gas 
in the nebula (-18 km s _1 ). The Hi structure is 10.0+1.7 pc in 
radius and expands at about 11 km s _1 . The associated atomic 
mass is 1500+500 M Q . 

Molecular gas with velocities in the range -27.2 and -14.0 
km s _1 surrounds the brightest parts of Gum 31. The sharp in- 
terface between the ionized and molecular material indicates 
that these gas components are interacting. We have estimated a 
molecular mass of (l.l+0.5)xl0 5 M . Two possible origins can 
be suggested for the molecular shell: 1) it could have originated 
in the collect and collapse process driven by the expansion of 
the H ii region; or 2) it could represent the remains of the natal 
cloud where the open cluster NGC 3324 formed. In both cases, 
a large amount of the neutral hydrogen in the H i shell may have 
originated in the photodissociation of the molecular gas. 

The volume density of the molecular cloud, higher than the 
electron density (33 cm -3 ), implies that the H n region is expand- 
ing. 

The emission in the far infrared correlates with the H ii re- 
gion and the H i envelope, indicating that the observed emission 
is probably related to the neutral and molecular envelopes. 



The distribution of the emission in the MSX band A, which 
closely delineates the H n region and correlates with the molec- 
ular emission, suggests the presence of a PDR at the interface 
between the ionized and molecular gas. 

A number of MSX, IRAS, and 2MASS point sources with 
IR spectra compatible with protostellar objects appear projected 
onto the molecular envelope, implying that stellar formation is 
active in the higher density cores of the molecular envelope 
around Gum 31, where massive stars or star clusters are prob- 
ably being formed. 

The optical image of the nebula does not show clear evidence 
of a central cavity, as expected in a stellar wind bubble, suggest- 
ing that the massive stars in the cluster have weak stellar winds 
or existed for a short period of time to develop an interstellar 
bubble in a high density interstellar medium. 
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